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Potential scattering problems governed by the time-dependent Gross-Pitaevskii equation are in- 
vestigated numerically for various values of coupling constants. The initial condition is assumed to 
have the Gaussian-type envelope, which differs from the soliton solution. The potential is chosen to 
be a box or well type. We estimate the dependences of reflectance and transmittance on the width of 
the potential and compare these results with those given by the stationary Schrodinger equation. We 
attribute the behaviors of these quantities to the limitation on the width of the nonlinear wave packet. 
The coupling constant and the width of the potential play an important role in the distribution of the 
waves appearing in the final state of scattering. 
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1. Introduction 

The nonlinear Schrodinger equation (NLSE), which has a cubic nonlinear term, 

i<f> t + <pxx + 2\<p\ 2 <p = 0, (1) 

appears in various fields of physics. ^ The NLSE can be derived as an amplitude equation of a system 
whose dispersion relation depends dominantly on the square of the amplitude. Among such systems 
is the envelope motion of coupled nonlinear oscillators with cubic interaction. 2) In nonlinear optics, 
both the self-focusing effect in two-dimensional (2D) systems and optical soliton propagation in one- 
dimensional (ID) systems are governed by the NLSE. 3 ' 4 ' Another important example is the Bose- 
Einstein condensed (BEC) system 5 ' 6) in which the macroscopic wave function of condensate atoms 
appears as the order parameter accompanied with the spontaneous breakdown of the U(l) gauge sym- 
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merry. 7 * In this case, the NLSE is regarded as the mean-field approximation of the Heisenberg equation 
for field operators and describes the time evolution of this macroscopic wave function with good ac- 
curacy. 5 ' 6 ^ 

One of the striking features of ID NLSE is its integrability. In particular, exact solutions un- 
der a given initial condition can be uniquely solved by the inverse scattering transformation (1ST) 
method. 8 ' 9) For the sech-type initial condition with suitable amplitude, the entire initial value problem 
is solved analytically and results in the /V-soliton solution. 10) However, the time evolution of the wave 
packet from the nonsoliton initial condition is relatively unclear since analytic expressions are rarely 
feasible. 

On the other hand, in the BEC system, it is natural to assume the existence of an external field to 
express the effect of gravity or quadratic traps for atoms. ' Thus, the term of external field is added to 
the conventional NLSE (1), and the equation is called the time-dependent Gross-Pitaevskii equation 
(TDGPE). So far, most analytic studies have assumed only linear or quadratic potentials, in which 
case the integrability of the systems are not spoiled and one can obtain analytical results by systematic 
application of the 1ST method. 12 ' In such cases, soliton initial conditions are also assumed. 

However, we can also consider spatially localized potentials where the integrability is manifestly 
violated. It is important to evaluate the role of nonlinearity on such potential scattering problems as 
the tunneling effect. In particular, time-dependent analysis of moving wave packets is intriguing since 
in situ observation of condensed atoms is possible in the BEC system, 5 - 6) although most studies deal 
with these kinds of problems as stationary ones. 14 ^ 

When we analyze the stationary potential scattering problems, we completely adopt the wavelike 
nature and the resonant phenomena brought about by the interference effect. On the the other hand, 
the spatially localized pulse is expected to exhibit the so-called wave packet effect in the scattering 
process. In addition, nonlinear effects are also interesting in the potential scattering problem. 

To clarify the definition, we shall use "the wave packet effect" and "the nonlinear effect" in the 
following senses throughout this paper, a) The wave packet effect means a phenomenon caused by 
the fact the wave packet is localized as a result of superposing many monochromatic waves, b) The 
nonlinear effect means a phenomenon that cannot be interpreted without considering the nonlinearity. 
To investigate the influence of these effects on the potential scattering problem, it is necessary to trace 
the dynamics of the system. Some authors have reported on this kind of problem assuming soliton 
initial conditions. 15 ' 1 * However, examples that take nonsoliton solutions as initial conditions have 
been rare because of extra complexities. In this paper, we numerically trace and examine the dynamics 
of the wave packets governed by ID TDGPE with the box- or well-type potential under the Gaussian- 
type initial conditions. 
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This paper is organized as follows. In the next section, the nonsoliton dynamics of the wave pack- 
ets without external field are analyzed. In §3, we evaluate and characterize the nonlinear wave packet 
in terms of the reflectance or transmittance, changing the magnitude of the nonlinearity, position of 
the initial wave packet, and the width of the potential. The section §4 is devoted to the discussion, and 
we interpret the results on the basis of wave packet and nonlinear effects. Extra complexities intrinsic 
to nonsoliton initial conditions are also argued in detail. The summary is given in the §5. 

2. Time Dependent Gross-Pitaevskii Equation and Scattering Problem 

In this section, we briefly summarize mathematical descriptions of the system to be considered. 
We restrict ourselves to the ID case throughout this paper. By virtue of scale transformation, we can 
put both coefficients <p t and <p xx of the TDGPE to be unity and we shall consider 

i<Pt = ~<f>xx + V(x)<f> + g\(f>\ 2 (f>, (2) 

where V(x) is an external potential applied to the system, and g is the coupling constant. To investigate 
the wave packet and nonlinear effects, the sign of g is important and we discard the possibility of 
g > throughout this paper. If we take g to be positive, which means a repulsively interacting field, 
the wave packet immediately expands and this rapid diffusing makes the amplitude of the wave packet 
very small. Therefore, the excitation of higher harmonic waves is extremely suppressed, and less 
manifestation of the nonlinear effect is expected. Moreover, these widespread wave packets share most 
of the scattering features with stationary plane waves in the linear limit. Therefore, we focus on the 
g < case, which means an attractively interacting field, in order to investigate the wave packet and 
nonlinear effects on the potential scattering problem. According to the theory of the partial differential 
equation, a finite and unique solution of eq. (2) exists for arbitrary initial conditions for the ID case, 
and instability or explosion of the solution observed in multidimensions never occurs even if we take 
g to be negative. 

Hereafter, we shall normalize the wave function <p as J^|0| 2 dx = 1. The initial condition of the 
wave packet is fixed to be the Gaussian type, 

<P(X,0) = -L e -I^) 2 + iv(x + , )^ (3) 

where -xq is the position of the center of the initial wave packet and v gives half of its velocity. We 
assume right-forward propagation of the wave packet, i.e., xq > and v > 0. The energy functional E, 
the Hamiltonian of the system, is defined as 

E= L ( l0x|2 + vixm2 + 2~ gl0 ' 4 ) dx - (4) 

Equation (2) can be derived straightforwardly from the Hamiltonian through standard canonical pro- 
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cedure. Since we have assumed g to be negative, E might take a negative value. In fact, for the initial 
wave packet located sufficiently far from the potential, the initial value of E becomes negative under 
the condition 

\ + v 2 + -£= < 0. (5) 

We consider the box- and well-type external potentials, 

V hox = 9{x) - 6{x - a), (6) 

Vwdi = -V o (0(*)-0(je-a)), (7) 

where a is the width of the potential and 9{x) denotes the step function. We define the reflectance and 
transmittance from these potentials as 



= lim f 

J-co 



R box = lim I \<f>\ 2 dx, (8) 



#weii = hm f \4>\ 2 dx, (9) 

*J—o 



= lim f 



\(f>\ 2 dx. (10) 

2+b 



The reason why we introduce b for the definitions of /? we ii un eq. (9) and T we \\ unh eq. (10) is as follow. 
For the well-type potential, part of the wave packet is trapped by the potential well and oscillates 
around the potential area. The distance of b is provided as a margin to distinguish the trapped portion 
and reflected or transmitted ones. The trapped portion never completely separates from the other parts 
of the wave packet, and continues to exchange a very small amount of their norms, and the limits 
in eqs. (9) and (10) do not exist in the strict meaning. However, for an evaluating measure, we use 
the values at t = 80 as if they were limiting ones. In the next sections, we vary g, xq, and a and 
investigate their influence on /?t>ox and T w& \\. For numerical integration, we employ the symplectic 
Fourier method 17 ' throughout this paper. 

3. Results 

In this section, the effects of nonlinearity on free propagation and potential scattering problems 
under the Gaussian initial condition are considered. 

3. 1 Free propagation 

In this subsection, we discuss the free propagation of a wave packet where no external potential 
exists. In this case, the 1-soliton solution of eq. (2) with V(x) = can be written as 

0(jc, t) = J—i]sech(xT] + 2ttf)e i{x t- t( i 2 -f 2)] , (11) 
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Fig. 1. Free propagating breather-like motion of a wave packet starting from the initial wave packet given in 
eq. (3) with x = 20, v = VT3, and g = -4. The nine wave packets show \(f>\ 2 at f=0, 2, . . . , 14, and 16 
from the left to the right. 



where 77 and £ are independent parameters and control the amplitude and velocity of the soliton, re- 
spectively. Once this form of solution is taken to be the initial condition, it never diffuses and mantains 
its own shape during time evolution. 

However, the Gaussian initial condition (3) leads to breather-like solutions. We show the time 
evolution of the wave profile \(f>\ 2 in Fig. 1. We also calculate the wave function in wave-number space 
(f> as 

1 r°° 

4>(k, t) = —— 0(jc, t)e lkx dx. (12) 
"V27r J-00 

We show a snapshot of |0| 2 taken at t = 16 in Fig. 2. In the wave-number space, the breathing motion 
is also observed and a notched structure grows on the surface of the wave packets. This structure 
is a result of repeated expansion and contraction in the wave-number space, i.e, expansion by the 
higher harmonic excitation and contraction by the dispersion effect (suppression of higher harmonic 
excitation). 

The breathing motion is a genuine nonlinear effect, since we would have observed that |0| 2 is 
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simply diffusing and |0| 2 remains still if g - 0. It is known that any solitary waves governed by eq. (2) 
with V(x) = finally split into a complete soliton part and a small oscillating tail (radiation) which 
rapidly leaves the soliton part in the long run. 18 ) Therefore, this breather-like behavior is considered to 
have a finite lifetime and the decay process is a rather transient phenomenon. Fortunately, this lifetime 
is sufficiently long to observe the breather-like motion. 

3.2 Box-type potential 

We move on to the main issue: the influence of the wave packet and nonlinear effects on the 
potential scattering problems. In this subsection, we consider a box-type (repulsive) potential (6). 

When we consider a stationary problem with g = 0, this is a textbook example of quantum me- 
chanics where the analytic expression for the reflectance is obtained. The expression reads 

4v 2 (v 2 - 1) 



1 + 



sin 2 (a Vv 2 - 1) 



(13) 



where v is twice the wave number of the incident plane wave. When sin(a Vv 2 - 1) = 0, R s becomes 
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Fig. 3. Reflectance from the box-type potential (6) for various values of g. The initial condition is the 
Gaussian-type wave packet (3) with xq = 5 and v = VT3, except for g — 0. The curve for g = corresponds 
to the linear case given by eq. (13). 




Fig. 4. Reflectance R^, m from the box-type potential (6) for various values of g. The initial condition is the 
Gaussian-type wave packet (3) with xq = 100 and v = VL5, except for g = 0. The curve for g = 
corresponds to the linear case given by eq. (13). 
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Fig. 5. Solid line shows |</>| 2 interacting with the box-type potential with width a = 0.5. The starting point is 
xq = 100 and g = -2. This is a snapshot taken at t = 26. The dashed line is for t = 0. 

and perfect transmission is realized. This is a kind of resonance scattering. 

Hereafter, we set x$ in eq. (3) to be 5 or 100 throughout this paper. The parameter v is also fixed 
at VT3. 

The dependences of 7?t>ox on g and a are shown in Figs. 3 and 4. The former is for jco = 5 and the 
latter, xq = 100. The curve for g = corresponds to the linear case given by eq. (13), the reflectance 
calculated from the stationery Schodinger equation. As mentioned above, the quantity R s becomes 
values as a increases. This is due to resonance and is also expected to occur periodically as the value 
of a grows larger. 

The behavior of R\, ox given by TDGPE (2) is drastically different. Firstly, maximum values of 
/?box for each g are totally suppressed for g < -2, although R\, ox is enhanced for the case of g = -2. 
Secondly, they never experience the perfect transmission resulting from the wave packet effect (This 
is a common feature of wave packet scattering. Such a result is expected even if we set g = 0, since the 
measure of the resonant monochromatic plane waves composing the wave packet is zero.), and they 
seem to be approaching their own constant values asymptotically accompanied with small oscillation 
as a increases, i.e., the periodic resonance structure is destroyed for g < -2 cases. 
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Fig. 6. Wave profile |</>| 2 colliding at the box-type potential. The parameters are the same as the ones used in 
the case ofFig. 5. 

This destruction should be regarded as the nonlinear effect because restoration phenomena of the 
resonance structure are expected if we use linear or weakly interacting wave packets starting from 
xo = 100. In fact, a wavy resonance structure seems to recover for the weakly interacting (g = -2) 
case after long free propagation, as shown in Fig. 4. This can be reasoned as follows: since relatively 
weak nonlinearity of g = -2 cannot prevent the wave packet from diffusing, it spreads and gains 
sufficient width for the plane wave approximation to be applied after long propagation. Therefore, the 
result approaches the linear case. 

To examine the validity of the plane wave approximation, we present Figs. 5 and 6. The former is 
the profile |0| 2 in the wave-number space interacting with the box-type potential with width a - 0.5. 
The starting point is xq = 100 and g = -2. This is a snapshot taken at t = 26. The latter is the wave 
profile \cp\ 2 in the real space colliding at the same box-type potential. The parameters are the same as 
the ones used in the case of Fig. 5. As we have shown in Fig. 5, |</>| 2 for g = -2 after propagation is 
concentrated around the center of the original profile. This means the spectral profile for g = -2 still 
keeps being localized. In contrast, the wave function for g = -2 in real space (Fig. 6) broadly spreads 
over the potential (a = 0.5). This implies the plane wave approximation is valid. Therefore, for the 
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Fig. 7. Typical wave shape including the portion trapped by attractive well-type potential (7) with a — 0.5. 
The wave packet located near the origin is the trapped portion. The initial condition is the Gaussian type 
wave packet (3) with xq — 5, v = VlT5 and g - -8. This figure shows the snapshot taken at t — 30. 
7\ ve u=0.233. 



linear or the sufficiently weakly interacting case, the behavior of /?t, ox for wave packets starting from 
xq - 100 is expected to be reminiscent of the plane waves and to exhibit the restoration phenomena. 

However, the restoration never occurs for strongly interacting wave packets (g < -2), which 
means the destruction of the resonance structure is due purely to the nonlinear effect. 

3.3 Well-type potential 

Next, we move on to the well-type (attractive) potential case. The potential is shown by eq. (7). 
The parameters are the same as those in the previous case, but the potential depth Vq is taken to be 10. 
For the stationary and linear case, the analytical expression for the transmittance is given as 



i + 100sin 2 (aVv 2 + 10)" 



(14) 



4v 2 (v 2 + 10) 

Perfect transmission is realized when sin(a Vv 2 + 10) = 0. However, in the nonlinear wave packet 
dynamics, we require more careful definition of the reflectance and the transmittance, since a certain 
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Fig. 8. Typical wave shape including the portion trapped by attractive well-type potential (7) with a — 5. The 
trapped portion forms a standing-wave-like structure in the potential well. The initial condition is Gaussian- 
type wave packet (3) with xq — 5, v — Vl-5, and g = -4. This figure shows the snapshot taken at t — 30. 
r well =0.640. 



amount of wave packets is trapped by attractive potentials. For example, snapshots of the trapped wave 
profiles are shown in Figs. 7 and 8. We can observe a standing- wave-like structure in the latter, and 
it swings back and forth in the potential area. Moreover, these trapping phenomena seem to be an 
intermediate state and the trapped parts continue to gradually emit a part of themselves mainly toward 
the left. Therefore, the reflectance defined by eq. (8) never converges, even after a very long time. 
Here, we employ expression (9) or (10) instead of eq. (8) to evaluate the wave packet and nonlinear 
effects. The value of margin b is chosen to be 30 in this work. 

The dependences of r we ii on g and a are shown in Figs. 9 and 10. The former is for xq = 5 and the 
latter, xq = 100. The curve for g = corresponds to the linear case given by eq. (14), the transmittance 
was calculated from the stationery Schodinger equation. The quantity T s becomes unity several times 
as a increases. This is also due to resonance and also is expected to occur periodically as the values of 
a grows larger. 

The main features resemble those of the box-type potential case. Firstly, the maximum values of 
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Fig. 9. Transmittance r we u over the well-type potential (7) for various values of g. The initial condition is the 
Gaussian-type wave packet (3) with x — 5 and v = VO. The curve for g = corresponds to the linear 
case given by eq. (14). 




Fig. 10. Transmittance T we n over the well-type potential (7) for various values of g. The initial condition is 
the Gaussian-type wave packet (3) with xo = 100 and v = Vl.5. The curve for g — corresponds to the 
linear case given by eq. (14). 
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Fig. 11. Portion trapped /V trap p e d by the well-type potential (7) for various values of g. The initial condition is 
the Gaussian-type wave packet (3) with xq — 5 and v = Vl-5. 



r we u for each g are totally suppressed for g < -2, although T we \\ is enhanced for the case of g - -2. 
Secondly, they never experience the perfect transmission resulting from the wave packet effect, and 
periodic resonance structure is destroyed for strongly self-focusing wave packets (g < -2). Thirdly, 
the wavy resonance structure seems to recover in the g = -2 case after long free propagation (Fig. 10). 
The reason for this restoration seems to be the same as that in the box-type potential case. 

Here, we evaluate the amount of trapped portions by subtracting the sum of reflectance (9) and 
transmittance (10) from unity, i.e., 

Xa+30 
I0| 2 dx. (15) 
30 

As mentioned before, these values are nothing more than estimates from the values at t - 80. Figures 
1 1 and 12 show the dependences of trapped on g and a. They basically show that the amount of trapped 
portion rises as \g\ and a increase, except for the relation between cases of g = -2 and -4 in the small 
a region. We attribute this phenomenon to a) stronger attractive interaction that contributes more to 
make the value of eq. (20) negative. In the next section, we explain that the negative value of eq. (20) is 
relevant to the trapping phenomenoa; b) wide potential width, where as it grows, the trapping capacity 
increases. Finally, we mention that the wavy structures seem to be accompanied with the resonances. 
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Fig. 12. Portion trapped/V tm pped by the well-type potential (7) for various values of g. The initial condition is 
the Gaussian-type wave packet (3) with xq = 100 and v = Vl.5. 



4. Discussions 

We studied the scattering problems of nonlinear wave packets as described in the previous section. 
One of the noteworthing properties of these problems is that the final reflectance or transmittance is 
a function of not only v but also the initial position of the wave packet. In §2, we showed strong 
modulation of the Fourier spectrum of a wave packet due to the nonlinearity which is shown in Fig. 2. 
The shape of the Fourier spectrum deforms and oscillates moment by moment during propagation. 
Therefore, the Fourier spectrum at the moment when the wave packet arrives at the potential area 
depends on the parameter xq, i.e., the distance between the starting position of the wave packet and 
the potential. The shape of the Fourier spectrum evidently affects the reflectance or the transmittance. 
This is the reason why the final results depend on x$. On the contrary, for the linear case, the initial 
Fourier spectrum is conserved under free propagation, and the role of xq is not important. 

The initial position of the wave packet also affects the result of the scattering problem through the 
alteration of the incident kinetic energy. The kinetic energy K is defined as 



f 



K= \<f> x \ 2 dx, 



(16) 



and the self-interaction energy as 



/ = \g f I0l 4 d*- (17) 
The breathing wave packet is always exchanging its kinetic and self-interaction energy even during 
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free propagation. As we can see from eq. (16), when the wave packet becomes steeper, the kinetic 
energy increases. Since total energy E = K + 1 is a conserved quantity in the case of free propagation, 
the negative self-interaction energy decreases to compensate the increase of kinetic energy. Therefore, 
the incident kinetic energy is also a function of the initial position of the wave packet xq. This incident 
kinetic energy directly fixes the wave number at the incident wave packet and becomes one of the most 
significant factors of the scattering problem. 

From the above considerations, any argument on potential scattering problems of a nonlinear wave 
packet requires the consideration of the initial position of the wave packet xq, except the soliton initial 
condition. In this work, we fixed xq to be 5 or 100. The reflectance and transmittance might be altered 
for a different choice of xq, while our main arguments are retained, i.e., the decay of the resonance 
structure and the existence of a trapped portion for well-type potentials might be observed. 

In the previous section, we showed the trapping effect for self-focusing wave packets by an attrac- 
tive potential. This phenomenon is interpreted as a purely nonlinear one. For linear quantum mechan- 
ics, this kind of phenomenon never occurs owing to the prohibition of energy level crossing between 
scattering and bound states. This is because time evolution of a wave packet from an initial state to 
a final one is fully described by a supeiposition of elements in the complete set of scattering state 
eigenf unctions, i.e., 



where A represents continuous energy eigenvalues and the integral is taken over all the scattering 
eigenstates. For nonlinear wave packets, however, the distribution of the conserved energy is always 
changing, as argued above, and a contribution from potential energy, 



Jo 

appears in the potential area. As seen from Figs. 7 and 8, the trapped parts of the wave packets are 
well separated from other parts, and hence it makes sense for us to define 



locally. If we calculate ea. (20), we actually get negative values. This trapped state can be considered 
a "dynamical bound state", since a part of the wave function is trapped and continues to oscillate 
dynamically, staying at around a potential having negative energy. 

It is worth mentioning that for the case of the box-type potentials, the reflectance of strongly self- 
focusing wave packets is observed to approach constant values for larger potential width a, and its 
dependence on parameter a disappears (Figs. 3 and 4). Because the norms of the self-focusing wave 
packets have finite values only in the narrow limited areas between potential ends and do not have 




(18) 




(19) 




(20) 
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sufficient extent to cover the whole potential area, the wave packets become insensitive to the opposite 
far end of the potential. In general, the resonance is a result of interference between the forward- 
propagating wave and the backward propagating one. However, the squeezed wave packets merely 
have a backward portion since they negligibly interact with the other far end of the potential. Such 
wave packets undergo the effective potential for large a, 

Vbox-eff = 0{X). (21) 

For well type potentials, the squeezed wave packets first fall off the cliff of the potential, 

VWeff-i - -V 9(x). (22) 

Then, they encounter the other side of the potential wall, and they effectively face 

Kvell-eff-2 = V 9(x - a). (23) 

The major part of the transmittance shown in Figs. 9 and 10 can be considered the remainders after we 
subtract the reflectance shown in Figs. 3 and 4 from unity. The reflected portion by the potential wall 
repeats reflection in the valley of the potential and is considered to constitute the dynamical bound 
states. 

5. Summary 

We numerically studied free propagation of wave packets governed by the TDGPE for various 
values of coupling constants g. The initial condition was taken to be the Gaussian form, which is 
different from the soliton solution. For the strongly self-interacting wave packets, diffusion in real 
space was suppressed and they exhibited breather-like behaviors. In wave-number space, the breathing 
motion was also observed, and a notched structure grew on the surface of the wave packet. 

We also numerically investigated the potential scattering problems under the same developing 
equation and initial conditions. The potential forms were chosen to be the box- or the well-type. We 
obtained the reflectance /?t>ox and the transmittance r we u for different values of coupling g and width 
of the potential a, and we compared them with the predictions made using stationary Schodinger equa- 
tions. We found that the reflectance R\, ox or the transmittance r we n is a function of the initial position of 
the wave packet xq. Other roles of nonlinearity are rather complicated, i.e., it sometimes enhances 7?t>ox 
or r we ii, but sometimes the opposite. However, there is a tendency that large \g\ decreases both Rb ox 
and r we n. For larger values of \g\ and a, R\, ox and T we \\ approach constant values and do not depend on 
a. 

We also observed the dynamically trapped portion of the wave packet. We estimated its amount, 
A^apped . by changing g and a and found that A^apped is an increasing function of g and a except in 
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the small g and a region. Whether this trapping effect is a perpetual or merely transitional one is not 
obvious and should be deteremined in future work. 

Finally, we make some remarks on the possibility of real experiments. The control of external 
environments is relatively easy in the BEC systems where we can confine condensate particles along a 
quasi-rectilinear line by tightening the laser beam trap. In addition, we can freely change the coupling 
constants by applying the Feshbach resonance technique. 19) Soliton-like pulses of BEC have already 
been created. 20 * If controllable local potential are realized, the possibility of observing and confirming 
our results by a real experiment is promising. 
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